Band-Limited Coronagraphs using a Halftone-dot process: 
design guidelines, manufacturing, and laboratory results 



P. Martinez 1 , C. Dorrer 2 , and M. Kasper 1 



ABSTRACT 

The Exo-Planet Imaging Camera and Spectrograph (EPICS) for the future 42-meter 
European-Extremely Large Telescope, will enable direct images, and spectra for both young 
and old Jupiter-mass planets in the infrared. To achieve the required contrast, several corona- 
graphic concepts - to remove starlight - are under investigation: conventional pupil apodiza- 
tion (CPA), apodized-pupil Lyot coronagraph (APLC), dual-zone coronagraph (DZC), four- 
quadrants phase mask (FQPM), multi-stages FQPM, annular groove phase mask (AGPM), high 
order optical vortex (OVC), and band-limited coronagraph (BLC). Recent experiment demon- 
strated the interest of an halftone-dot process - namely microdots technique - to generate the 
adequate transmission profile of pupil apodizers for CPA, APLC, and DZC concepts. Here, we 
examine the use of this technique to produce band-limited focal plane masks, and present guide- 
lines for the design. Additionally, we present the first near-IR laboratory results with BLCs that 
confirm the microdots approach as a suitable technique for ground-based observations. 



Subject headings: instrumentation: high angular resolution 
tion 



techniques: high angular resolu- 



1. Introduction 



B y the end of 2018, challenging projects such as EPICS (Exo-Planet Imaging Camera and Spectro- 
graph, Kasper et al. J 20081)') for the future 42-meter European-Extremely Large Telescope, or PFI (Planet 
Formation Imager, llvlacintosh et alJl2006bl) for the Thirty Meter Telescope (TMT) will enable direct images, 
and spectra for warm self-luminous and reflected-light Jovian planets. These instruments will operate with 
eXtreme Adaptive Optics sy stem (X AO) designed for high S trehl ratios (i.e. ~ 90% in //-band), as for 
SPHERE feeuzit et al.ll2008l) and GPI (IMacintosh et a l. 2006a), forthcoming planet-finder instruments with 
first light planned in 201 1. 

Several coronagraph concepts have been studied extensively in the past years, with the objective of 
finding optimized desi gns that can s ufficiently suppress the on-axis starlight, allowing f ai nt companions 



direct detection (e.g. iMalbetl (119961) : ISivaramakrishnan et all (1200 ll) : iGuyon et all (120061) : Martinez et al. 
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(I2008ah ). Among them, the band-limited coronagraph (IKuchner & Traubll2002l . BLC) has been proposed 
to completely remove starlight. The BLC has the advantage o f being less sensitive to the primary m irror 
segmentation, unavoidable with ELTs, than for other concepts ([Si varamak rishnan & Yaitskoval 120051) . and 
to provide achromatic behavior. Additionally, in [Martinez et al.1 (12008 ah . we pointed out the interest of 
such concept on an ELT for either very bright object detection, or for the search of planets at large angular 
separations. 



Several BLC proto t ypes have been developed during the past y ears (e.g. lDebes et aU2004UTrauger et al. 
2004l : lcrepp et allboodJTrauger & Traubll2007l : [Moody et all2008l) for visible wavelength application. Sev- 



eral technical approaches have been used: (1/) gray-scale pattern written with an high-energy beam sensi- 
tive glass (HEBS) using e-beam lithography, (2/) Notch filter pattern - binary mask - written with thick 
Chromium layer on a substrate, dry-etched with high density decoupled plasma, (3/) gray-scale pattern 
manufactured with vacuum deposited metals and dielectrics. Even if electron-sensitized HEBS glass can 
accurately accommodate continuous range of transmission, the darkening of the HEBS glass under electron 
bombardment is accompanied by a determined phase shift, while the technique suffers from a lack of expe- 
rience in the near-IR. Same constraints apply for the vacuum deposited metal technique. The notch filter has 
the advantage of being intrinsically achromatic. These designs, consisting of a particular implementation 
of small structures (stripes of opaque material with width of about some microns) must be finely c ontrolled 
in size, spacing, and opacity. Mask errors and tolerance are discussed in lKuchner & Spergell (120031) . where 
requirements might be strong for near-IR application. 

In this paper, we examine the use of a halftone-dot process, namely mic rodot technique, to repro duce 
a continuous mask profile as already done for pupil apo dizer for SP HERE (|Martinez et al.ll2009allbh . and 
being manufactured for the JWST NIRCam coronagraph (lKristll2009l) . These masks consist of distributions 
of opaque square pixels (called dots) to reproduce the continuous transmission of a filter with several advan- 
tages: relative ease of manufacture, achromaticity, reproducibility, and ability to generate continuous trans- 
mission ranges, withou t introducing wavefront errors. Besides, mask errors can be easily pre-compensated 
(borrer & Zuegelll2007l) . 



Section 2 describes the microdot design principle and properties, while Sect. 3 provides guidelines 
for the design. Section 4 presents monochromatic, and polychromatic results obtained in laboratory in the 
near-infrared. Finally Sect. 5 concludes on the suitability of the microdots approach for producing BLCs in 
the context of ground-based instruments. 



2. Principle 



In its general scheme, a microdots filter is an array of dots (i.e. pixels) that are either opaque or transpar- 
ent. It is fabricated by lithography of a light-blocking metal layer deposited on a transparent glass substrate. 
To be st approximate specificati ons, the dot distribution is regulated by a Floyd-Steinberg dithering algo- 
rithm (IFloyd & Steinberg! 1 1976b . based on error diffusion principle. This allows the accurate gen eration of 
grey levels, and rapidly varying shaping functions. This algorithm is presented, and discussed in Juiichney 
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19871 Il988h for printing technique applications, as well as for laser beam shaping (IDorrer & Zuegell l2007). 



Additionally, this algorit hm was successfully u sed for producing microdots pupil apodizers for the apodized 



pupil Lyot coronagraph ( Martinez et al. 2009ailb[) . 



By principle, the microdot mask - the BLC amplitude function - is sampled, subject to the dot size 
(see Fig. \Q. This issue will be further discussed in Sect. [3] from a design specification point of view. In 
the following, we will examine its impact on the coronagraphic effect by comparing to that of a continuous 
transmission mask. 



We consider a teles cope with aperture function A, and a band-limited amplitude mask function from 
Kuchner & Traubl J2002h : 

M(r) = N(l -sinc(— (1) 



where A is the wavelength of the application, r the radial coordinates in the image plane, D the telescope 
primary diameter, e the bandwidth which rules the inner-working angle of the coronagraph (IWA hereafter), 
and finally N is a constant of normalization insuring that M(r) ^ 1. 

The Fourier transform of M{r) denoted ^{u), is therefore: 

^M-w(«(£)-3*n(£)) (2) 

where u states for the radial coordinate of the pupil plane, 5 the Dirac-function, and IT the top-hat function. 

In a BLC, the incident electromagnetic field propagates from the telescope aperture to a focal plane where 
the mask, M(r), is applied. After the occulting mask (filtering for the low-frequencies), the pupil is re- 
imaged in a second pupil plane where a pupil-stop is placed (high-frequency filter). Then the scientific 
image is recorded in a second focal plane where the detector is installed. The coronagraphic effect can 
be well described in the second pupil plane (i/f pp , hereafter), as the convolution product of ^{u) by the 
telescope aperture A(u), times the pupil-stop function denoted S(u), such as: 

\pp P - \_Ji(u) ★ A(u)] x S{u) (3) 

where the star (★) states for the convolution product. The power of band-limited coronagraph comes from 
the properties that ^(u) is equal to zero everywhere but \u\ < Sj, i.e. the power spectrum of M(r) has power 
in a limited domain of frequencies. As a result, in the second pupil plane (Eq. O, the convolution product 
of ^#(w) by the pupil aperture confines the diffracted light in the vinicity of the pupil edges as exhibits in 
Fig. [2] (left, where the pink color states for numerical noise level). This energy can be completely removed 
with the pupil-stop, S (w), s imilar to the aperture fun ction, but reduced in diameter. A development of the 



BLC theory can be found in lKuchner & Trau b (2002), although we note that the present formalism assumes 
that the amplitude mask is purely real, i.e. it does not take into account the relationship between opacity and 
phase shift in real physical material, which impacts the performances. 

A microdot BLC will have a different power spectrum than a continuous mask because its function is sam- 
pled by the dots, which forces its Fourier transform to be periodic, and the algorithm used to distribute dots 
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introduce high-frequency noise. As a result, the convolution product of ^ (u) by the pupil aperture still 
confines the diffracted light in the vinicity of the pupil edges, but some energy remain inside and outside the 
geometrical pupil as exhibits in Fig. |2](right). The residual part inside the geometrical pupil, not filtered by 
the pupil-stop, will therefore set a limit in the achievable performance. The intensity of the residual energy 
as seen in Fig. |2] (right) is 4 to 5 order of magnitude bellow the maximum intensity of the ring surrounding 
the pupil aperture. The amount of residual energy in the pupil plane is a function of the dots size, and the 
bandwidth of the function (e), since these two parameters directly impact the sampling of the design. This 
point will be further analyzed in Sect. I3.2l with simulations. 

Although BLC designed with microdots can not achieve perfect cancellation of starlight, it is important 
to note that a ground-based instrument does not require ideal performance of a BLC, but can accomodate 
with relaxed performance since limitation will first come from the quasi-static speckle halo level. Planet- 
finder instruments will use speckle calibration strategies to improve contrast delivered after correction of the 
turbulence by an Adaptive Optic system (AO). If a simultaneous speckle rejection techniqu e is used (e.g . 
simu ltaneous differential im aging, spectral deconvolution...), as selected for SPHERE (|Beuzit et al.l l2008). 
GPI (Macintosh et al.l l2006ah . and EPICS, performance of a coronagraph does not necessarily need to be 
better than this quasi-static speckle halo level (10~ 5 to 10~ 6 ). A deeper contrast is then achievable through 
appropriate data reduction. 



3. Mask design guidelines 

3.1. Specified functions 

We consider two BLC masks with identical function (Eq. |2]), but different bandwidth (i.e. different 
e values, see Fig. [3]>. Hereafter, BL5 and BL10 correspond respectively to e - 0.17 (IWA = 5A/D) and 
e = 0.086 (IWA = 10A/D). Even if BL10 has a large IWA (~ 0.08 arcsec in H-band for a 42-meter 
telescope), only suitable to the search of companions at large angular separ ations, it has a smoother transition 
in between the low and high-transmission parts of its profile, and can be compared with BL5 in the light of 
the technique employed (i.e. sampling issue). 



3.2. Size of the dots 

The selection of the dots size is in first approximation - as discussed in Sect. [2]- formally equivalent 
to a sampling problem. Better results are obtained with smaller pixels, since this allows finer control of the 
local transmission. This is particularly an important issue in region of the mask where transmission is very 
low (center of the mask). Since BLC masks are placed in the focal plane of an instrument, their dimensions 
are usually expressed in A/D units. Therefore, a simple metric to quantify the sampling of the function is the 
ratio of the Airy unit by the dot size (p). Optimally, this metrics must be defined at the shortest wavelength 
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for which the mask is designed to operate, and can be expressed as: 



F# x A m j n 

s = (4) 



where F# is the /-number on the mask and A m i„ the shortest wavelength of the application (i.e. insuring 
that s remains greater than the specification for all wavelengths). Using simulation maps of microdots BLC 
(assuming specified dot spatial distribution), we analyzed how the dot size affects the coronagraphic perfor- 
mance with respect to continuous idealistic masks. Our simulations make the use of Fraunhofer propagators 
between pupil and image planes, which is implemented as fast Fourier transforms (FFTs) generated with 
an IDL code. Dots are sampled by one pixels to prevent from excessive size of array critical in the high s 
values domain. Validity of the dot sampling has been verified with simulations assuming different sampling 



(1, 4 and 16 pixels per dot), as discussed in Martinez et al . (2009a). An effect is only expected in the very 



low s domain at high-spatial frequencies, out of the field of view of interest. These monochromatic, and 
aberration-free simulations were addressed for the two band-limited configurations we manufactured (BL5 
and BL10). Results are presented in Fig. |4l Note that a continuous mask provides a total on-axis source 
suppression (hence not plotted in the figure) assuming the idealistic conditions of these simulations. As 
expected (Sect. |2]), ideal performance of BLC can not be met (i.e. contrast deeper than 10~ 12 , in Fig. H]). 
However, several s configurations (8, 16, and 32) yield interesting contrast levels, deep enough to meet con- 
trast requirements of planet-finder projects. Considering our application (H-band, A$= 1.6 microns, with a 
contrast requirements of 10~ 6 at 5A/D after the coronagraph, e.g. EPICS), suitable configurations are s - 8 
and s = 16, as the s = 32 situation yields to very small dots size (2.5 microns) comparable to the size of 
the wavelength of light. In this particular situation, it is difficult to predict how the field will react to such 
grating, where scalar diffraction theory does not apply. Additionally, since e impacts the gradient from the 
very low-transmission maintained at the center to the high-transmission part (bandwidth of the function), 
better performance are obtained when e is larger (Fig. |4j left and right). The issue is that the very low 
transmission maintained at the center of the BLC functions is significantly impacted by the dot size. This 
point will be further discussed in Sect. |~~ 



3.3. Function bandwidth 

The bandwidth parameter of BL functions (e) directly impacts the gradient from the low to the high- 
transmission part of the profile, and rule the IWA of the coronagraph, it is therefore important to analyze its 
effect on the performance. We simulated several cases with different e values (from 0.095 to 0.8, i.e. IWA 
from 1 to 9 A/D, see Fig. [5]> for the 5 = 16 case. Results gathered in Fig. [5]show that for IWA larger than 
3A/D (e < 0.28) performance in the halo are not £-dependent, while only a slight difference is observable 
on the peak. For IWA in between 3 and 2A/D, although an important degradation of performance appears, 
coronagraphic capabilities are not negligible (contrast of 10~ 5 at IWA). Under 2A/D IWA, situation degrades 
further. 
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3.4. Dot opacity 

Whatever the regime of wavelength for which the device is designed, the optical density (OD) is a 
critical issue to avoid leaks of the starlight in the center of the detector image. The OD is guaranteed by 
sufficient thickness of the metal layer (denoted e hereafter). For visible application, Chrome layer are com- 
monly used, while for near-IR applications, more opaque material are mandatory (Aluminum, for instance). 
The OD is defined at a specific wavelength, and for a given material, and thickness such as: 

0« = ^ (5) 
A x ln(W) 

where k describes the linear attenuation of the optical wave (material-dependent), e stands for the propaga- 
tion distance (i.e. the material thickness), and A is the operating wavelength. Equation [5] is further detailed 
in appendix |AJ The OD must be therefore carefully defined, accordingly with expected performance of the 
coronagraph (s and e-dependent, Figf4]and[5]). 



4. Experiment 
4.1. Optical setup 

The optical setup is a near-infrared coronagraphic transmission bench developed at ESO (Fig. [6]). All 
the optics are set on a table with air suspension in a dark room and are fully covered with protection panels 
forming a nearly closed box. The entrance aperture is full-filled, and has a 3 mm diameter (<1>) made in a 
laser-cut, stainless-steel sheet to an accuracy of 0.002 mm. BLC masks were installed at an F/48.4 beam. 
Re-imaging optics were made with IR achromatic doublets. The pupil-stop has a diameter D stop - 0.84 x <J>, 
and remains the same during the experiment. Its optimization gave emphasis to the achromatic behavior 
of BLCs. The quality of the collimation in pupil planes was checked and adjusted using an HASO 64 
Shack-Hartmann sensor. A pupil-imager system was implemented to align the pupil-stop mask with the 
entrance-pupil mask to complete alignment in the x- and y- directions. The coronagraphic focal plane was 
localized using a visible CCD mini-camera with a HeNe laser light and tuned in the final IR image on the 
detector. The IR camera used (the Infrared Test Camera) uses a HAWAII Ik x Ik detector, cooled to 105 K 
degree with a vacuum pressure of 10~ 5 mbar. Internal optics were designed to reach a pixel scale of 5.3 mas 
(almost 8 pixels per A/D). Experiment was done in H-band using either a narrow filter (AA/A = 1.4%), or a 
broadband filter (AA/A = 24%). The Strehl ratio of the bench was evaluated to ~ 92%. It was determined 
by measuring the peaks intensity ratio of the experimental PSF to that of the theoretical PSF normalized 
to the total intensity. The theoretical PSF is created through two different methods, both converging to 
the same Strehl ratio: (If) by performing the foward fast Fourier transform (FFT) of the autocorrelation 
of an oversampled and uniformly illuminated entrance pupil image from our telescope pupil mask, (2/) by 
performing the FFT of a simulated aperture function with radius determined from the experimental PSF on 
the basis of photometric criteria. 
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4.2. Prototypes 

Four prototypes were manufactured by Precision Optical imaging (Rochester, New-York), two BL5, 
and two BL10 masks. All the masks were designed for 1.64/mi (H-band), and fabricated using wet-etch 
contact lithography of an Aluminum layer {OD = 8+, e = 2000A) deposited on a BK7 substrate (A/ 10 
peak-to-valley, 0.5 inch diameter). Antireflection coating has been applied for each faces (R < 0.5% from 
1.2 to 1.8/mi). For each design (BL5, and BL10), two dot sizes have been produced: 5 and 10 microns (i.e., 
s = 16, and s = 8) denoted by (a) and (b) respectively (see Fig. [7]). The OD has been specified to guarantee 
against leakage greater than the intrinsic limitation of the microdot technique (for s = 16, and s = 8, see Fig. 
|4|. All the masks have been inspected and cleaned up before integration (Fig. [7]). Profiles were measured 
at 1.0/rni. The spatially resolved transmission has been obtained after background subtraction and flat field 
normalization. Profile accuracy is of about 5% of the specification. The error is mostly localized in the 
outer part of the mask (high-transmission part), while the center part (for the low-transmissions) is highly 
accurate (Fig. [8]). The error in the outer part took its origin in a calibration issue of the process, that will be 
corrected for with new prototypes. 

4.3. Results 

The following metrics are used to evaluate coronagraphic abilities to suppress the on-axis starlight: 

• t stands for the total rejection rate: ratio of total intensity of the direct image to that of the corona- 
graphic image. 

• To stands for the peak attenuation rate: ratio of the maximum of the direct image to that of the 
coronagraphic image. 

• 'lo stands for the contrast: ratio of the coronagraphic image at a given angular separation to that of the 
maximum of the direct image to the intensity, azimuthally averaged. 

4.3.1. Data acquisition and reduction 

To achieve high-dynamic range measurements, a serie of 3s short exposure images averaged over 3 
mn, and the presence or absence of neutral density filters, were employed. Neutral density filters were only 
applied on non-coronagraphic images. PSFs and coronagraphic images were dark subtracted by turning off 
the light. Data reduction process correct for bad pixels, background, and scales images by the exposure time 
and optical density. 
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4.3.2. Results with BL5 

BL5 prototypes have been tested with the narrow (AA/A = 1.4%), and broadband (AA/A = 24%) filters. 
Coronagraphic images are presented in Fig. |9j while profiles of BL5(a) and BL5(b) are gathered in Fig. \\0\ 
(top). All coronagraphic runs made with BL5(a) and BL5(b), whatever the filter, roughly yields to identical 
performance. BL5 demonstrate an achromatic behavior as expected for such coronagraph. No difference in 
between BL5(a) and BL5(b) has been observed, therefore regarding to the performance obtained, s = 16 and 
s = 8 are suitable configurations. A discrepancy of ~1 order of magnitude for both BL5(a) and BL5(b) has 
been revealed on the peak, while in the halo it evolves between 2 and 3 order of magnitude (see theoretical 
results as a function of s in Fig. HI left). Despite t his discrepancy these f irst results (gathered in Table [J) are 
better than the SPHERE prototypes performance (Boccaletti et al.ll2008h . and experimental results reached 



on the same bench with an APLC ([Martinez et al.ll20 09a). However, we note that these prototypes (Four 



Quadrant Phase Mask, and APLC) were tuned for very small IWA (1, and 2.3 A/D respectively compared 
to 5 and 10 for the BLC described here). Contrast evolves from ~ 3 x 10~ 5 at IWA to ~ 3 x 10~ 8 at 20A/D 
(cut-off frequency of SPHERE-like AO system). 



4.3.3. Results with BL10 

Results presented with BL10 prototypes (Fig. [101 bottom) correspond to AA/A = 24%. As for BL5, 
we did not notice any difference between performance reached either by BL10(a), or BL10(b). Contrast 
evolves from ~ 1 X 10~ 7 at IWA to ~ 1 x 1(T 8 at 20A/D, while the peak attenuation rate is about 10 5 . A 
discrepancy of about one order of magnitude has been revealed on the peak as for BL5, while in the halo 
it evolves between 1 and 2 order of magnitude at IWA, and at 25A/D depending on which configuration of 
BL10 we are looking at (Fig. |4] left). As expected BL10 results went beyond BL5 performance, although 
the IWA is quite large. The coronagraphic image exhibits speckle halo (Fig. [9]) induced by the non-perfect 
nature of optics used in the bench. This point is further discussed in the following section (Sect. 14.4.2b - 



4.4. Performance limitations 

The discrepancy revealed in the experiment may find its origin in several error sources, among them - 
for the most important ones -, the scatter-light due to imperfect optical components that create speckles in 
the scientific image, and mask profile errors. Note that the source diameter (0.082/l/D), and the alignment 
errors of the masks are negligible considering the IWA of the BLCs, while the impact of the filter bandpass 
is not considered because of the achromatic behavior of BLCs (confirmed in the experiment). We also 
neglect the fact that the pupil-stop is not black-coated (e.g. black anodization) which may produce potential 
reflections. 
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4.4.1. Mask profile errors 

Simulations using real measured profiles (i.e. global profile measured after free-space propagation 
of the BLC masks and using an imaging optical system, Fig. [8] for instance), showed that profile errors 
are actually responsible of one order of magnitude loss on the peak and while the halo is less affected. 
Assuming their respective measured profile errors, the rejection factor for BL5 and BL10 (s = 16) are 
respectively limited to: r mJ t = 6400, and T mas k = 108000. 

Ideally, going to very small pixel size improves the accuracy of the mask profile transmission (i.e. the 
sampling problem). In practice, two constraints set a limit on the interest of a very small pixel size: (1/) 
getting good accuracy becomes more difficult because fabrication errors become more important as the pixel 
size decreases , (2/) when the pixel size is comparable to the wavelength of light, the transmission is affected 



by plasmons ([Genet & Ebbesenl 120071 : iHuang & Zhu 2007). In such situation, the transmission might be 



strongly dependent on the wavelength. Therefore, it is very likely that efforts must be focused on reducing 
fabrication errors to avoid an increase of transmission, as a result of a reduction of the metal dots during 
the isotropic wet-ech lithography process. An improvement of the pre-compensation algorithm applied on 
the digital design, correcting for edge effects on the light-blocking metal dots, by estimating the feature size 
that would be obtained after fabrication is thus critical. 



4.4.2. Wavefront errors 

The wavefront error (wavefront high-frequency components) of the optical components prior to the 
pupil-stop imposes an important limitation when reaching high-contrast. The optical components on the 
bench upstream of the pupil-stop have a standard quality: flats mirrors of A/10 ptv (i.e A/35 rms), and 
achromatic doublets of A/5 rms, over the full diameter (25.4 mm) at 633 nm, which corresponds respectively 
to ~ A/300 and ~ A/ 45 rms for the 3 mm pupil diameter (assuming the quality scales linearly with the beam 
size). Taking into account all components, and assuming that phase added quadratically, it corresponds 
to an overall approximated total amount of wavefront error of ~ A/61 rms at 1.6 microns, (i.e. 24 nm 
rms). With simulations using a Fraunhofer propagator IDL code, and assuming theoretical mask profile, 



and u sing typical power-law of manufactured optics, described by PSDs with / variation (Du parre et al. 



2002), where / is the spatial frequency, we actually found that 12 nm rms prior to the pupil-stop, already 
imposes a contrast floor in between 10~ 7 and 10~ 8 , which is in good agreement with the wavefront error 
estimation discu ssed above, which is commonly used as a rule of thumb for order of magnitude estimation, 



as performed in iRiaud et all (|2003h . For 12 nm rms, the total rejection factor for BL5 and BL10 (s = 16) 



are respectively limited to: T wave f ront = 40400, and T wave f ront = 188200. 
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4.4.3. Performance estimation vs. experimental results 

Assuming that these independent errors are added quadratically, the global rejection factor of BLs including 
mask error, and wavefront errors is: 




which correspond to tbls = 6300 and tblio - 93600. These estimations are in agreement with experimental 
results (Table [j}. Note that x refer to a global performance estimation (integrated metric of the overall 
images, see Sect. 14. 3b . and can therefore not inform for a local coronagraphic profile behavior. Figure [TT] 
gathers experimental and theoretical results for BL10 (s=16) with a broadband filter (AA/A = 24%), with 
simulated profiles obtained assuming independently mask profile error and wavefront error (12 nm rms), as 
well as a profile including both of these error sources in the same time. Even if it is difficult to ascertain 
source of errors, the profile error is consistent with the discrepancy found at small angular separation (i.e. 
peak level error), while the wavefront errors is probably the dominant source of error in the halo. In that 
situation, the intrinsic limitation of the microdots technique is therefore not a limiting factor. 



5. Conclusion 

We have described the development and laboratory experiment of Band-limited coronagraphs using 
a microdots design in the near-infrared. In this paper, we provide design guidelines, and demonstrate the 
microdots technique as a promising solution for BLC for ground-based observations. 

We have shown with numerical simulations that although total starlight cancellation is not possible, 
theoretical contrast offer with the microdot approach are deep enough to guarantee that the BLC will not set 
a limit on the performance of a ground-based instrument. We note that the theoretical treatment presented 
in this study does not consider the complexity introduced by potential spectral, and polarization effects of 
the physical mask. 

We identified two sampling configurations suitable for near-IR experiment: s = 16 and s - 8, yielding 
to identical performance in the experiment. Additionally, we pointed out that the interest of the microdots 
technique in the light of contrast and IWA requirements is not dependent of the function bandwidth (e) 
assuming TWA > 3A/D. This already meets EPICS (20-30 mas in H-band), or SPHERE (0.1 arcsec in H- 
band) standard requirements. However, we note that stronger requirements of the IWA (1 or 2 A/D, ultimate 
goal of SPHERE), will set a limit on the interest of the technique. 

With prototypes we have demonstrated impressive performance, where limitations have been presum- 
ably (on the basis of simulations) introduced either by mask error, or wavefront error of the bench. Improve- 
ment of the results presented and resumed in Table Q] are foreseen - at least for the peak attenuation - with 
a new set of prototypes that will provide a better accuracy of the profile in the outer part of the function 
(calibration issue of the manufacturing process). 
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Additionally, these raw polychromatic results presented belong as the first tests of BLC in the near- 
IR, and were managed without active wavefront correction, nor particular data reduction post-processing. 
Performance reached in the experiment went beyond to the SPHERE prototype performances, APLC results 
obtained so far with a microdo t apodizer on the same bench, and yield to similar contrast than the ones 
presented in ICrepp et al.1 (|2006) (using 4 th order notch-filter mask, in a monochromatic visible-wavelength 
domain) although IWA are not exactly identical. 

Ultimately, these final prototypes will be implemented on HOT (the High-Order Testbench, the AO- 
facility at ESQ, |A ller Carpenti er et al.l 120081) . and being compared with others (FQPM, APLC, and Lyots, 
Martinez et al.ll2008bl l2009ah w ith atmospheric turbu lence generator and AO-correction, as already initi- 
ated with intensive simulations ( Martinez et al.ll2008ah . Results of this experiment will be presented in a 
forthcoming paper. 

Although this study was carried out in the context of R&D activities for EPICS, it is potentially appli- 
cable to upcoming instruments such as SPHERE, or GPI. We note that the interest of the technique presented 
in the paper for space-based operations, is subject to science cases (IWA and contrast requirements). For 
instance, high star-planet contrast (10~ 10 in the visible) at less than 0.1" (Terrestrial planet), will very likely 
required a wavefront stability only reachable with a space telescope. In that situation, the technique em- 
ployed here will impose a limit on the accessible contrast. Notch-filter masks might be more appropriate. 
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A. Optical density of a surface layer 

A lossy dispersive medium can be described for optical propagation by a complex index n + ik where n 
and k are related to the linear susceptibility. The index n describes the evolution of the phase of the optical 
wave and k describes the linear attenuation of the optical wave. The propagation of an electromagnetic wave 
of wavelength A = 2nc is described by a complex wavevector ^ (n(A) + ik(A)), where both the index and 
attenuation are wavelength-dependent. 

The input and output electric fields after a propagation distance (e) are related by: 

E ou t = E in x exp L^(n(A) + ik(A)e\ (Al) 

The intensity transmission (T - E out /Ei n ) for a thickness e and wavelength (A) is therefore: 

/ Ank{A)e \ 

T - exp (A2) 



The optical density is by definition the base- 10 logarithm of the inverse of the intensity transmission, i.e. the 
natural logarithm of the intensity transmission divided by the natural logarithm of 10: 

and is therefore dependent on the layer thickness (e), the operating wavelength (A), and the factor k (material- 
dependent). 
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Fig. 1. — Simulation maps of a circular band-limited mask. Top: Scan of the center (grey level pattern). 
Bottom: corresponding optimized microdots pattern. Grey-scale color evolves from 0% (black) to 100% 
(white) transmission. 




Fig. 2. — Example of pupil plane images before effect of the pupil-stop, with a continuous BLC (left), and 
a microdot BLC (right). The arbitrary false color distribution has been chosen to enhance contrast for the 
sake of clarity. The scale evolves from noise-level (black / pink colors) to maximum intensity (red-color). 
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Metrics 


T 


to 




^12X1 D 


^20^/D 


^25A./D 


BL5 (a) 
BL5 (b) 


2410 
1804 


2554 
2038 


3.7 x 10~ 5 
3.0 x 10" 5 


5.6 x 10~ 7 
8.1 x 10" 7 


2.7 x 10~ 8 
2.6 x 10" 7 


5.3 X lO" 8 

8.5 x 10~ 8 


BL10 (a) 
BLIO (b) 


81606 
83177 


97375 
131809 


1.5 x 10~ 7 
1.3 x 10" 7 


7.7 x 10" 8 
6.3 x 10" 8 


3.7 x 10~ 8 
4.3 x 10" 8 


1.3 x 10~ 8 
2.1 x 10~ 8 



Table 1: Summary of coronagraphic results obtained with BL5 and BLIO (AA/A - 24%). 




20 40 60 

Angular separation in X/D 



Fig. 3. — Band4imited functions used in simulation, and in the experiment. 
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Fig. 4.— 



Impact of the dot size on the coronagraphic performance for BL5 (top) and BL10 (bottom). 
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Fig. 5. — Impact of the mask bandwidth on the IWA (left) and performance (right). 
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Fig. 6. — Schematic setup of the coronagraphic test-bench. 
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Fig. 7. — The four prototypes manufactured as seen in their integration mounts (right), successive Shadow- 
graph inspections (middle and top-left, x 50 and x 100), as well as a spatially resolved transmission of BL5 
(b) recorded at 1.0 microns (bottom-left). 
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Fig. 8. — Measured profiles of BL5 (top), BL10 (bottom), and specification for the two dot-size prototypes 
(5 and 10-micron). 
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Fig. 9. — From left to right: PSF image (with a ghost, right side of the image), BL5 coronagraphic image, 
and BL10 coronagraphic image (AA/A = 24%). The image dynamic has been chosen to enhance contrast 
for the sake of clarity. 
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Fig. 10. — Experimental results obtained with BL5 (top) and BL10 (bottom) prototypes, all profiles are 
azyimuthally averaged. 
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BL10 (s=16 & AX/A = 24%) 
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Fig. 11. — Resume of BL10 (s=16 and AA/A = 24%) theoretical (green curve) and experimental results 
(red curve) compared to simulation profiles obtained, when including error sources such as mask profile 
error (dotted black curve), wavefront error (12 nm rms, full-filled black curve), or both of these errors (blue 
curve). 



